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A technique to deal with Coulomb electron distortions in the analysis of 
(e,e'p) reactions is presented. Thereby, no approximations are made. The 
suggested technique relies on a partial-wave expansion of the electron wave 
functions and a multipole decomposition of the electron and nuclear current 

•i— I . 

rS | in momentum space. In that way, we succeed in keeping the computational 



times within reasonable limits. This theoretical framework is used to calculate 

the quasielastic (e,e'p) reduced cross sections for proton knockout from the 

valence shells in 16 0, 40 Ca, 90 Zr and 208 Pb. The final-state interaction of the 

ejected proton with the residual nucleus is treated within an optical potential 

model. The role of electron distortion on the extracted spectroscopic factors 

is discussed. 
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I. INTRODUCTION 

For a long time it has been recognized that the exclusive (e,e'N) reaction in the quasielas- 
tic (QE) region is a powerful tool for studying the single-particle motion inside the nucleus, 
and is a testing ground for the different available nuclear models. One of the principal in- 
terests in the exclusive (e,e'N) reaction is to extract the nucleon spectral function P(p, E) 
from the cross section. This spectral function can be interpreted as the joint probability to 
remove a nucleon with momentum p from the target nucleus and to find the residual system 
at an excitation energy E. Related to these spectral functions, spectroscopic factors and 
occupation numbers are often studied. They are a measure for the validity of the indepen- 
dent particle model (IPM). The spectroscopic factor S n ij m (E) gives the probability to reach 
the single-particle state specified by the quantumnumbers nljm in the residual nucleus at 
an excitation energy E. The occupation number N n ij m gives the number of nucleons in the 
single-particle state nljm in the target nucleus and involves an integration of the spectro- 
scopic factors over the complete excitation energy range Jl|. In the IPM the states above 
(under) the Fermi level are completely empty (filled) and the total hole (particle) strength 
is situated at a fixed single-particle energy. The deviation from full (no) occupancy for the 
orbits below (above) the Fermi level is a measure for correlations neglected in this mean-field 
approach. 

The occupation probabilities in even-even nuclei have been calculated within several the- 
oretical frameworks. Most models go beyond the mean-field approach and partially account 
for short- and/or long-range nucleon-nucleon correlations Jy-H]. Occupation probabilities 
for the single-particle states which considerably deviate from the IPM value were obtained. 
Moreover, it is demonstrated that the single-particle hole strength is fragmented over a broad 
range of energy. In particular, occupation numbers for the proton 3sl/2 orbit in 208 Pb have 
been calculated varying from 1.42 to 1.66 pointing towards a strong depletion of this 
hole state in the ground state of 208 Pb. From an experimental point of view, the CERES 
method H was developed in an attempt to obtain absolute occupation numbers from experi- 



mental data. The model uses only relative spectroscopic factors and allows to account, in an 
approximate way, for the strengths at high missing energies, not accessible for experiment. 
With this method, the 3sl/2 occupation number in 208 Pb is found to be 1.57(10). 

Although the advantages of the quasielastic (e,e'N) process to study spectroscopic fac- 
tors are widely recognized, the extraction of these factors from experiment is still not free of 
ambiguities. For example, depending on the model used in the analysis of the 208 Pb(e,e'p) 
reaction, the spectroscopic factor for the transition to the groundstate in 207 T1 ( 3sl/2 
hole) varies from 0.40 || to 0.71 fK|. A reliable determination of spectroscopic factors 
requires an accurate knowledge of the (e,e'N) reaction mechanism (photoabsorption mecha- 
nism, final-state interaction (FSI) of the ejected nucleon with the residual nucleus) and the 
exact treatment of the Coulomb distortion of the scattered electrons, especially for heavy 
nuclei. 

In this paper we present results from systematic calculations of (e,e'p) cross sections for 
a number of even-even target nuclei and various kinematical conditions and confront them 
with data taken at NIKHEF. The extracted spectroscopic factors are compared with the 
corresponding values deduced within other theoretical approaches pH|^r5|. Much attention 
is paid to the effect of electron distortion on the calculated cross section. It is pointed out 
that, especially for scattering off heavy nuclei, an exact treatment of these effects is highly 
needed in order to reproduce the shape of the measured cross sections and, consequently, to 
obtain reliable spectroscopic factors. 

This paper is organized as follows. In section II the theoretical formalism for the (e,e'N) 
reaction is outlined. The derivation of the cross section is divided in two subsections treat- 
ing the electron and the nuclear aspect of the (e,e'N) reaction. The technical details are 
dealt with in appendix |A]. The numerical details of the adopted approach are discussed in 
section III. The formalism is applied to electro-induced one-proton knockout reactions from 
a number of medium-heavy target nuclei in section IV. Finally, some conclusions are drawn 
in section V. 



II. FORMALISM 



A. Cross section 



In this paper we describe the process in which an electron with four-momentum k(e, k) 
and spin polarization m Sk is scattered from a target nucleus at rest with a rest mass Ma- 
The detected electron is characterized by its four-momentum k'(e', k') and spin polarization 
m s i . The energy transfer to the nucleus u = e — e' is supposed to be sufficient to eject a 
nucleon N (proton or neutron) with four-momentum pn{En,pn) and spin projection m SN 
out of the target nucleus leaving the residual nucleus with four-momentum pb{Eb,Pb)- The 
differential cross section and the Feynman amplitude rrifi for this process are related as 
follows: 



dV 1 



e' 2 \p N \E N J2\m fl \ 2 S(u - S N - E x - E N - E B + M N + M B ) . (1) 
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Throughout this paper we adopt natural and unrationalized Gaussian (a = e 2 ) units. In this 
relation Sn stands for the separation energy of a nucleon out of the target nucleus and E x 
denotes the excitation energy of the residual nucleus. The rest masses of the ejected nucleon 
and the residual nucleus are given by Mn and Mb- The angles Q e (9 e ,(p e ) and Qn(@n,(Pn) 
specify the scattered electron and ejected nucleon with respect to the chosen reference frame. 



At this point this reference frame is not further specified. The sum J2ij implies a summation 
over all final states (electron and nuclear) and an average over the initial states (electron and 
nuclear). We only have to sum over these final states which satisfy the energy conservation 
relation. 

In the Born approximation the transition amplitude rrifi can be written in terms of 
matrixelements of the electron J^J and nuclear J^ ucl charge- current four- vector in momentum 
space in the following way: 

1 /■ 1 

m fi = ~2^2 J d ^2_ \-\2 + i S < /el J el,/.(-g)l^ e >< fn\ J"nucl(<T ) K > ■ (2) 



The initial and final electron states are denoted by \i e > and |/ e >. The target nucleus and 
final nuclear state consisting of a residual nucleus and an ejected nucleon are represented by 
\i n > and \f n >. 

The Feynman amplitude ra/j can further be rewritten as follows 

m / 1 = ^~2 / d <l{ 7^12 < fe\Pel(~q )|*e > < /n|Pnucl(<?) Kn > 




+ 



^ (-1) A « < / e |Jel,A,(-? )l»« >< fn\Jnncl,-X q (q)\in > 

\q=±l 



(3) 



c<j 2 — | q | 2 + irj 

The spherical components of the electron and nuclear current operators are taken with re- 
spect to the rotating reference frame (x q , y q , z q ) (Fig. p] (a)). In this way the third component 
of the current operator is directly related to the charge operator through the charge-current 
conservation relation. 

B. The leptonic part 

In this section we elaborate on the electron matrixelement < fe\J&,n(—Q)\ie > i n the 
expression for the Feynman amplitude. The relativistic electron charge-current operator in 
coordinate space reads 

J°(f) = _ e vj) e t(r) 4f*(f) 
J A (f) = -e \i> et (r) a * 

with \l/ e (r) the electron field operator in coordinate space. The initial and final electron 
wave functions are defined according to 

<f\i e >=^l(f), (5) 

<f|/ e >=*!(r), (6) 

and stand for four- dimensional Dirac spinors. They are solutions of the stationary electron 
Dirac equation: 

(a.(-iV) + (3m e + V) fcj^ (f) = e ^ (f) , (7) 



(4) 
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where m e is the rest mass of the electron and V is the scattering potential. The additional 
quantumnumber m Sk uniquely determines the electron wave function. 

Dealing with high-energetic electrons the electron mass can be neglected with respect 
to its total energy and the Dirac equation can be written down in the ultrarelativistic limit 
(e = \k\). In the Dirac- Pauli representation for the a and (3 matrices and in the absence of 
an external potential V the solutions of equation (|7|) are given by : 






(f) =u e (k,m Sk ) e 



ik:\ 
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V2 
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x\ll W, 



vf<(^ 



ik.f 



(8) 
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The spinors xll (^fc) can be expressed in terms of the Pauli spinors and the matrixelements 



of the Wigner Da -matrix, i.e., 



xlil^k) = £x^>)^m s >feA,o) . 



(9) 



The angles Qk = (9k, ipk) specify the momentum k with respect to the chosen reference 
frame (x,y,z). The Wigner T>2(R k ) matrix represents the rotation of the reference frame 
(x, y, z) over the Euler angles Rk = (<fk, 8k, 0) in the basis spanned by the eigenvectors of 
the operators S 2 and S z . 

Assuming a central potential V = V(r), the electron wave functions are evaluated by 
a phase shift analysis based on a partial-wave expansion. Indeed, the Dirac Hamiltonian 
(H = a.k + V(r)) commutes with the angular momentum operators J 2 and J z and with the 
operator K = f3 {a.L + 1} but not with the orbital momentum operator L 2 . As such, we 
derived a complete set of operators with common eigenf unctions represented by \&^ m (r): 



J 2 Kjmif) 



e \§> e . If) 



Jz KiJr) = m * 



(10) 
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Kjm\ 
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We can construct the partial waves ^ e K j m (r) as follows: 
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Kjm(r) = %m(r) 



with 



/ ^yk"!^ r ,a) 



V 



^y{^ r ,a. 



\ 

) 



(11) 



1 = 3-1/2 if « = -(; + 1/2) 

/ = j + 1/2 if k = j + 1/2 . 



We introduce the common notation I 



l=3 + \ 



1 = 3- 






(12) 



The spherical spin-orbit eigenspinor yfy 2 (p, r , a) is defined in the following way 
y n n / 2 (Q r ,a)= £ <lm l l/2m s \ ] m>Y lmi (tt r )x 1 J l 2 M) ■ 



(13) 



mim s 



Each partial wave (11) can be easily proved to satisfy the eigenvalue equations (pi]) under 



the condition that the radial electron wave functions G\Ar) and FfAr) are solutions of the 
following second-order differential equations: 



&Gh(r) + %®&4=GfAr) 



dr 2 ^lj^ 
+ 



E-V{r) dr ^lj y 



(E - V{r)f - am + fi ™fl GL(r) = 



(14) 



dr 2 ^Ijv) + E-Y(r) dr J '^ 



r £-V(r) 

dV»/dr d cw >> 

dr X ' ' ■ ' ' 

For each partial wave Ij the second-order differential equation for G\Ar) has to be solved 
numerically. For the regular solutions one imposes the following boundary conditions: 



+ 



F£(r) = 



limG£(r) = 0, 

i — >o J 

Km — GL(r) = for />0, 
r->o dr J 

and one obtains the corresponding solution for Ff.{r) through the relation: 

GUr) = (l-J)FUr). 



(15) 



(16) 



The asymptotic behaviour of the radial electron wave functions for Coulomb potential 
scattering are given by (k = |A;|) 



, r e ( s (1 A sin(fcr - Zvr/2 + ^ (tot) - r?ln(2A:r)) 

hm Gy(r) = (/-/) 1 , (17) 

sinffcr - Ivr/2 + <£: e (tot) - n ln(2Ar)) 
lim F£(r) = 1 { . (18) 

The phase shift 5^ e ° reflects the influence of the scattering potential V. It consists of two 
parts, i.e., the Coulomb phase shift of and an additional phase shift S^ € . For a Coulomb 
potential generated by the Z protons in the nucleus, the Coulomb phase shift is defined 
according to (rj = —Ze 2 ): 

<rt = aigr(! + i + in). (19) 

Due to the fact that the scattering potential V is spin- independent, one can easily verify 
that the total phase shift is /-independent, i.e., 

^e,e (tot) _ ^e,e (tot) _ ^c,e (tot) s^Q) 

Finally, the electron wave function \p£ (r) is expanded in terms of these partial waves 

Ijm 

The initial and final electron wave functions have to satisfy the outgoing (+) respectively 
incoming (— ) boundary conditions. Knowing the asymptotic behaviour of the radial electron 
wave functions, the coefficients a l ^ are fixed by 

C fe(±) = E vi sm3k (R k )^ l e ± ^ it0t \l-l)Y* mi (Q k ) < l mi l/2m s \jm > . (22) 

m 3 mi \ *• 

At this point only the scattering potential V remains to be specified. In general the 
central Coulomb scattering potential generated by Z protons is given by 

\ j-r poo 

V(r) = -AuZa- / p(r')rdr' - 4nZa / p(r')r'dr' , (23) 

T JO Jr 

with p(r) the nuclear charge density normalized according to 4tt / °° p(r)r 2 dr = 1. In the 
forthcoming discussion we have taken this charge density to correspond with a homogeneous 
spherical charge distribution of Z protons within the nuclear radius R. 
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By switching off the scattering potential V one can easily verify that the solution 



coincides with the free electron wave function (|8|) since the differential equations (0) reduce 
to the differential equations for the spherical Bessel functions. In this way a sensitive testing 
case for our numerical approach is found. 

We want to stress that the problem of Coulomb distortion of the initial and final electron 
in the electron scattering process is solved to all orders. Earlier work in this field by Boffi 



et al. [|TIJ handled the electron distortion in an approximate way through a high-energy 



expansion of the electron wave functions combined with an expansion in powers of Za. 



The DWEEPY code (JIT] used in the analysis of the NIKHEF data adopts this approximate 
treatment of electron distortion. To lowest order in Za it was proved that electron distortion 
effects could be approximated by an effective momentum approach (EMA). This means that 
the plane wave in eq. (|J) has to be replaced by 



ik.r ik .r 



e — * ~r e > ( 24 ) 

K 



with 



:cff . 3Za, 



k ={k + ~2R )Sk ■ (25) 

Clearly this approach is very easy to handle and worth comparing with the complete distorted 
wave approach so that its degree of accuracy can be estimated. 

C. The nuclear part 

In a previous paper H], we have shown that at low values of the missing momentum, 
meson-exchange currents (MEC) and long-range effects only slightly affect the calculated 
(e,e'p) cross section. As we will restrict ourselves to QE (e,e'p) reactions at low missing 
momenta only the one-body part of the nuclear four-current is retained. Hereby we adopt 
the operator as dictated in the non-relativistic impulse approximation: 

Pnuci(r)= Y, eGW^Sir-n), (26) 

i=l...A 



i=l...A 



■4d(r)= £ E ~M (W~ ri) + 8(r- r t )V„ 



This nuclear charge-current four-vector refers to A non-interacting point-like nucleons with 
mass Mj. To correct for the finite extent of the nucleons, the Sachs electromagnetic form- 
factors Ge and Gm are introduced. 

As for the electron wave functions, the final nuclear wave function is determined through 
a phase shift analysis after an expansion in partial waves. The final nuclear state is taken to 
be a linear combination of one particle-one hole excitations \C; ujJM > out of the A-particle 
groundstate \i n > with C = {h, p}. The hole state h is characterized by the quantumnumbers 
nh,lh,jh an d energy e h . The continuum particle state is specified by the quantumnumbers 
p = (l,j) and the energy e p = E^ — M^. The isospin nature of the particle-hole state is 
denoted by t g . The particle-hole state in the coupled scheme is defined according to 



\C;luJM> = J2 <] h -m h jm\JM>{-l)^- m ^\ph-\uo) > , (27) 

with the uncoupled particle-hole state defined as follows 

\ph- l (u) >= c+(e p )c h \i n > , (28) 

and u = e p — €h- The operators c + and c denote single-particle creation and annihilation op- 
erators. The radial wave functions for the bound hole states are solutions of the Schrodinger 
equation with a Hartree-Fock potential generated with an effective interaction of the Skyrme 
type (SkE2) |15|. The continuum particle states are evaluated within an optical potential 
model (OPM) [pJ]] . The physical radial wave functions are regular in the origin and behave 
asymptotically (r — > oo) according to 

i i \ r-*op f 2^ s ' m ( k p r - l ' 7T / 2 -V ln 2fc p r+<5 ; "' £:p(tot) ) 
Ppv) ► y Kkj, r € P > U ' 

(29) 
Mr) ^ e h < 
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where 77 and the momentum k p = \k p \ stand for 

kp = 2fi N e p with /j, N = M N (A — 1)/A the reduced mass of the nucleon , 

(oUj 

_ (Z-l)afj, N 
'I— k p 

The complex phase shifts caused by the nuclear and Coulomb part of the optical potential 
are denoted by 5?f p and of (8?f p{tot) = Sff p + of). 

Given the asymptotic behaviour for the radial single-particle wave functions and imposing 
that the ejected nucleon wave function satisfies the incoming boundary conditions, the final 
nuclear state \f n > is given by 



\ fn> = Y. H 4vr ^ \ T < ihni h jm I JM >< lmi-m s | jm > 

Ijmmi JM V l ^ K V l 

. r ll,€p(tot) 

xe-^ Y* m[ (n N )\(l h3h ,l 3 y,LuJM> . (31) 

In order to derive this expression the target nucleus is considered to be a spherical nucleus 
in the J n = + groundstate. In addition, the residual nucleus is described by a pure hole 
state h with respect to this target nucleus groundstate. 



D. The Feynman amplitude 

As the initial and final electron wave function and the final nuclear state are expanded in 
partial waves, it is common to decompose the electron and nuclear charge-current operators 
in the Coulomb, Electric and Magnetic multipole operators of rank JM (q = \q\): 

Tf M (q) = - fdr V x [jj(qr)yy {Jtl) (a r )} J(r) , 
q j 

TZ S (q) = Jdr 3j (qr)yy {JA) (a^J(r) , 
Mf M \q) = Jdrjj(qr)Yj M (n r )p(r) , (32) 

with the vector spherical harmonics defined according to 

yf{L,i)(tt) = E < LM L l\\ JM > Y LML (n)e x , (33) 

\M L 
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and e\(\ = 0, ±1) the standard spherical unit vectors corresponding with the unit vectors 
(e x , e y , e z ) in the (x, y, z) reference frame (Fig. [J] (a)). 

Accordingly, in momentum space the charge- current operators can be written as: 

p(q) = ^J2i J Y* JM (n q )M^(q), 

JM 

J Xq (q) = -yfa £ t J J [TjiM + WZ S (q)} V J MX {R q ) (34) 

J>1,M 



with J = \J2J + 1 and the Euler angles R q = (<p q , 9 q , —4> q ) defined in Fig. [I] (a). 

It is well-known that when neglecting electron distortion effects, the differential (e,e'N) 
cross section can be written in terms of four structure functions containing all the nuclear in- 
formation. In such a Distorted Wave Born Approximation (DWBA) approach each structure 
function is multiplied with an analytical factor containing the leptonic information. This is 
no longer valid in the Coulomb Distorted Wave Born Approximation (CDWBA) approach 
as the electron part can no longer be separated from the nuclear part. Consequently, when 
accounting for Coulomb distortion effects one has to perform a multipole expansion for both 
the electron and nuclear charge-current operators. 

Combining equations (§) and ([34]) and applying some basic properties of the Wigner 
T> J (R q ) matrices the Feynman amplitude m/, reads as (q^q 11 = uj 2 — q 2 ) 

™fi = E (-^tSi/ 00 ^ [< feWruti^e >< /«|M£SL(?)|i» > 



(2tt) 3 Jo 

Q {< fe\T e L ^(q)\i e >< fn\Tl^ L mn 



LM L 

-,2 



q^ + 17] 

+ < f e \r^ L {q)\i e X f n \T^ ML {q)% >}] . (35) 

The superscript e and n refer to the electron and the nuclear multipole operators. We have 
deliberately chosen to work out the leptonic and nuclear matrixelements in momentum space. 



Earlier electron distortion calculations by Jin et al. [|I(J and Udias et al. [12 1 evaluate the 
transition matrixelements in coordinate space. In order to make their calculations feasible 
the nucleon formfactors are evaluated at the asymptotic value q — k — k! . The major 
advantage of our approach is that the momentum-dependence of the nucleon formfactors 
can be handled exactly. 
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Here we will solely calculate the unpolarized (e,e'N) cross section ([!]) so we need to 
evaluate 

EK*I 2 = ^ E E K*l 2 - (36) 

ij Z m Sfc m Sfc/ m B m SN 

The summation over the initial and final states involves a sum over the initial and final 
electron polarizations and a sum over the polarizations of the recoiling nucleus and the 
ejected nucleon. In appendix |A] this Feynman amplitude is further worked out. Summarizing 
from appendix [A], one can state that the calculation of the (e,e'N) cross section is reduced 
to the evaluation of a large number of leptonic radial integrals 7^Lj U - 2 (e, e'; q) and a set of 
reduced transition matrixelements C(C; qui J) containing all nuclear information. We stress 
that the technique developed here can be easily extended to polarization processes. 

III. NUMERICAL PROCEDURE 

In order to derive the exclusive (e,e'p) cross section (||) in the CDWBA we need to 
evaluate the Feynman amplitude rrifi discussed in the previous section and appendix [XJ. The 
numerical procedure is schematically sketched in Fig. |2|. From a numerical point of view 
the evaluation of this transition amplitude is cumbersome as it involves an integration over 
the complete q range and two infinite sums, i.e., the sum over the different multipolarities 
L in the multipole expansion of the leptonic and hadronic current and the sum over the 
angular momentum j% originating from the partial-wave expansion of the scattered electron 
state. Angular momentum selection rules make sure that the other summations in the 
equations (|A3|) and (|A10|) have a finite range for fixed values of j\ and L. 

When accounting for electron distortion effects, the integrandum in the integral over 
q peaks at the effective momentum transfer g eff = \k cS — k' cS \. As the EMA is only an 
approximation of electron distortion effects the integrandum is spread around this value and 
the integration in g-space has to be performed in an interval [g m i n , <?max] around q eS . It is 
worth noting that in the absence of electron distortion effects, the integral over q vanishes 
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and the standard DWBA expressions are retained. The integrandum then reduces to a 
(^-function representing the momentum conservation relation q = \k — k'\. 

The finite extent of the nucleus puts a constraint on the number of multipolarities L which 
have to be retained in the multipole expansion of the nuclear current given in eq. fl34"|). In 
the calculations we systematically observe convergence when including multipolarities up to 
I max ~ 2qR where R denotes the radius of the considered target nucleus. In a similar way the 
number of electron partial waves which contribute to the (e,e'p) cross section, is restricted 
by an upper limit ji, m ax- It can be easily verified that the electron partial waves GfJr) and 
Fu (r) corresponding with large values for Ij are negligible for values of r within the nucleus 
range. For that reason, these electrons can cause no nuclear transitions. Consequently, to 
a required accuracy, only a finite number of the electron partial waves contributes to the 
electron scattering cross section. The number of electron partial waves actually contributing 
to the cross section depends also on the electron energy. The higher the electron energy the 
more partial waves will be required. The numerical evaluation of the (e,e'p) cross sections is 
getting complicated due to the large number of electron partial waves to consider. This is a 
result of the long-range character of the Coulomb interaction. The limit Z — ► (equivalent 
with turning off the electron distortions) can be considered as a severe test of the accuracy 
of the numerical techniques and a convergence test for the electron partial waves. For Z = 0, 
the electron wave functions reduce to plane waves. Accordingly, the DWBA cross section 
should be retained. As will be demonstrated in the forthcoming sections, our code has been 
checked to comply with this requirement. 

Another important feature of our CDWBA approach is that the radial integrals TZlj 1 j 2 
(|A4|), which are the heart of our numerical procedure, do not depend on the scattering angles 
6 e and 9 P . Consequently, our numerical procedure is optimized for calculating the (e,e'p) 
cross section for these specific kinematical conditions where the electron and proton scat- 
tering angles are varied and the other electron characteristics are kept fixed. The complete 
missing momentum range of the (e,e'p) cross section for proton knockout from the different 
hole states can then be calculated with a stored set of radial integrals. 
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IV. RESULTS 

Up to now, most of the high-resolution (e,e'p) experiments performed at NIKHEF (Am- 
sterdam), Saclay, Mainz and MIT-Bates have been carried out by using either parallel or 
constant q — uo kinematics. Both correspond with in-plane experiments: the ejected pro- 
ton is detected in the scattering plane spanned by the initial and final electron. In parallel 
kinematics the proton is detected in the direction of the momentum transfer. By varying 
the incoming e and outgoing e' electron energies or/and the scattering angle 9 e , different 
values for the momentum transfer q — k — k! and consequently the missing momentum 
Pm = Pp — (fare reached. For constant q — u kinematics the energy-momentum transfer is 
kept fixed and the proton angular distribution is measured. The missing momentum is de- 
fined positive when the ejected proton lies in the half-plane of the initial electron momentum 
and bordered by the momentum transfer. In the other half-plane the missing momentum is 
negative. 

Most of the experimental data are presented in terms of the reduced cross section ex- 
tracted from the measured cross section in the following way (p m = \p m \, p p = \p p \) 

1 dV 

Pm\Pmi Ar) = ^ , , , n , n , ,-, (37) 

PphpUep de'di l e a\ l p ah p 

with <j ep the off-shell electron-proton cross section. We stress that only in the plane-wave 
impulse approximation (PWIA) the reduced cross section coincides with the nucleon spectral 
function P(p m ,E x ), i.e., the probability to eject a nucleon with momentum p m from the 
target nucleus while leaving the residual nucleus at an excitation energy E x . As soon as the 
FSI, electron distortion and many-body nuclear currents effects come into play this quantity 
can no longer be interpreted as the nucleon spectral function. In comparing our (e,e'p) 
results with the available data we have divided the calculated cross sections with the ccl 



prescription |16j for a ep . The same procedure was applied to the experimental cross sections 
presented in this paper. Moreover, the calculated curves are scaled with a spectroscopic 
factor which accounts for the fragmentation of the single-particle strength. 
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The results of our model calculation are compared with the predictions from three other 
model calculations. Firstly, we confront our results with the non-relativistic CDWBA model 
of Bom et al. DTTJ . This model is at the basis of the DWEEPY code often used in the analysis 
of the NIKHEF (e,e'p) data. In the latter model the FSI is treated in a non-relativistic 
optical potential calculation similar to ours. In contrast with our model, the bound state 
wave functions are calculated in a Wood-Saxon well. The rms radius of the bound state 
wave function is fitted to reproduce the shape of the measured reduced cross section and 
the well depth is adjusted to reproduce the experimentally observed separation energy. In 
our calculation, we use the bound state wave functions as obtained from a Hartree-Fock 
calculation with a density- dependent effective interaction. Accordingly, in our approach the 
spectroscopic factor is the only parameter adjusted to the data. Concerning the treatment 
of electron distortion effects the two models are very different. Whereas in our calculation 
Coulomb electron distortion effects are treated to all orders, the CDWBA model of Boffi et al. 
implements electron distortion effects within the high-energy expansion as briefly mentioned 
in the theoretical discussion of section |11B| . In comparing the results obtained with these two 



non-relativistic models one can study in how far an exact treatment of electron distortion 
effects is required in the analysis of (e,e'p) reactions. 

Our results for the reduced cross sections and the corresponding spectroscopic factors are 
also confronted with the completely relativistic calculations of Jin et al. [jTOj and Udias et 



al. |T2[ . In line with our approach, the two models handle the electron distortion in an exact 
distorted wave calculation. The main difference with our model occurs in the description of 
the photoabsorption process and the initial and final nuclear system in the (e,e'p) process. 
Jin et al. and Udias et al. work in a totally relativistic framework. The bound state wave 
functions are calculated from the Dirac equation with a scalar and vector potential which 
are parametrized fits to relativistic Hartree potentials. The wave function of the knocked 
out nucleon is the solution of the Dirac equation with a relativistic optical potential. Two 
different prescriptions for the relativistic nuclear current operator are considered. They are 
referred to as the ccl and cc2 current operators and follow the conventions of ref. [[16] . 
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A. Parallel kinematics 

In this section we deal with the quasielastic (e,e'p) reaction from 16 0, 40 Ca, 90 Zr and 208 Pb 
in parallel kinematics. In Table | we specify the studied kinematical conditions. They all 
correspond with measurements performed at the NIKHEF electron accelerator. The (e,e'p) 
cross section for these different target nuclei are calculated in the CDWBA framework as 
outlined in the previous sections. The FSI of the ejected proton with the residual nucleus is 
handled within an OPM. For the medium-heavy nuclei 40 Ca, 90 Zr and 208 Pb the potential 
as derived from the Schwandt parameterization |2!| is considered. This optical potential is 



known to provide a good description of the elastic (p,p') scattering data over a large range of 
target mass and incident proton energies. The target nucleus 16 is out of the range of nuclei 
used in the parameterization of this global optical potential. Therefore, for the 16 0(e,e'p) 
calculations, we adopt the optical potential which is directly extracted from a recent analysis 
of elastic 16 0(p,p') scattering data at T p = 100 MeV and use the parameterization quoted as 
"WS" in ref. [[17]]. In order to study the effect of Coulomb distortions the (e,e'p) predictions 
from the DWBA and CDWBA model are compared. We stress that these two models only 
differ in the way the Coulomb distortions are described. In the DWBA they are completely 
neglected, whereas in the CDWBA they are treated exactly. 

The CDWBA reduced cross sections for electro-induced proton knockout from the lpl/2 
and lp3/2 shell in 16 are confronted with the NIKHEF data in Fig. |. The DWBA and 
CDWBA curves for each state are multiplied with one and the same spectroscopic factor. 
This spectroscopic factor is extracted from a least-square fit of the CDWBA reduced cross 
section to the data. The multiplication factors as extracted from our calculation adopting 
the WS optical potential are given in Table ||. Table |1| also lists the spectroscopic factors 
obtained within the non-relativistic CDWBA model of the Pavia group |TT[ as reported 



in ref. [17|. Comparing the results presented in this work and the predictions outlined in 



ref. [17], a similar degree of agreement with the data is reached. The extracted spectroscopic 



factors agree within 10%. From Fig. R] it is clearly seen that the calculated reduced cross 
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sections reproduce the measurements very well and electron distortion effects, although 
small, improve the agreement with the data especially for knockout from the lp3/2 orbit in 
16 0. 

We also performed calculations for electro-induced one-proton knockout from the ld3/2 
and 2sl/2 shell in 40 Ca. The results are plotted in Fig. [|. Electron distortion effects seem 
to follow the same pattern as observed for electron scattering from 16 0, but the effect is 
now more pronounced. From the 16 0(e,e'p) and 40 Ca(e,e'p) results one can already trace 
the main effects of electron distortion on the reduced cross section in parallel kinematics: 

• Electron distortion shifts the reduced cross section towards higher missing momenta. 
This can be explained by considering that a virtual photon exchanged between the 
electron and the nucleus will carry a momentum q eS instead of q (q eS > q). From 
equation (p5|) and the definition of p m one deduces that this shift will be decreasing 
with increasing p m . 

• The shape of the reduced cross section is mainly modified at the minima and maxima. 
Clearly, electron distortion not only manifests itself in an effective momentum shift 
but also in a focusing effect of the electron beam onto the target nucleus. 

The curves in Fig. ^ are scaled with a spectroscopic factor obtained from a least-square 
fit of the CDWBA results to the data. In Table [If J we compare the spectroscopic factors 
from our analysis with those obtained from the non-relativistic analysis with the DWEEPY 
code |ll| and those extracted from the two complete relativistic calculations by Jin et al. 
and Udias et al. @. 

Comparing the spectroscopic factors for proton knockout from the ld3/2 shell two main 
features can be observed: 

• The spectroscopic factors obtained with the non-relativistic models are in very good 
agreement with each other but are considerably smaller than the relativistic values 
obtained with the cc2 nuclear current operator; 



• The spectroscopic factors extracted within the relativistic models seem to be very 
sensitive to the prescription for the off-shell nuclear current operator. The cc2 current 
operator results in a spectroscopic factor for the 2d3/2 state that differs with more 
than 20% from the ccl result. The ccl current operator is obtained from the cc2 
current operator using the Gordon decomposition and should produce similar results 
for on-shell nucleons. 

The appreciable difference between the relativistic and non-relativistic approaches is rather 
surprising considering that the proton kinetic energies dealt with are typically of the order 



of 100 MeV. According to Jin et al. ]23J and Udias et al. fT2|j24}| the noticeable difference 
between the relativistic and non-relativistic spectroscopic factors is caused by the stronger 
absorptive part in the relativistic potentials. Even though all optical potentials reproduce 
the elastic proton-nucleus scattering data to a more or less similar degree, the quenching of 
the reduced cross section due to the final-state interaction of the ejected proton with the 
residual nucleus can differ by 15% adopting a relativistic or non-relativistic optical potential. 
This can be attributed to the behaviour of the optical potential in the nuclear interior. One 
could however doubt whether the interior part of the optical potential can be constrained 
in elastic proton scattering processes that are typical surface events. 



Hedayati-Poor et al. [25| attribute the difference between the relativistic and non- 
relativistic spectroscopic factors to the nuclear current operator. They show that the 
non-relativistic reduction of the relativistic transition amplitude results in an effective non- 



relativistic current operator which depends on the strong scalar and vector potentials |26j for 
the bound and the continuum single-particle states. Instead of using this medium-modified 
non-relativistic nuclear current, we adopt the standard non-relativistic nuclear current op- 
erator in our calculations. In our opinion, this is justified as long as the sensitivity of the 
relativistic results to the choice of the relativistic nuclear current operator is not cleared up. 
Concerning the spectroscopic factors obtained for proton knockout from the 2sl/2 shell 
(see Table [111]), the different models give very different predictions. In conformity with 
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the calculation of Udias et al. we describe rather poorly the reduced cross section around 
p m = 0, especially in the negative missing momentum region. This results in a spectroscopic 
factor which is not very reliable. However it has to be stressed that, in contrast with what 
was done in the analysis of ref. ||18|| , no attempt has been made to improve the results by 
adjusting the parameters of the optical potential and/or by adjusting the bound-state wave 
characteristics (rms radii and binding energies). 

The next target nucleus we considered is 90 Zr. The calculations cover knockout from 
the different valence shells in 90 Zr for two different proton kinetic energies (T p = 70 and 
100 MeV). We investigate to what extent the reduced cross sections for knockout from the 
outermost shells (2pl/2, 2p3/2, lf5/2) are affected by electron Coulomb distortion effects. 
Secondly the results of the complete calculation are confronted with the available data. In 
Fig. [| the reduced cross sections derived within the DWBA (neglecting electron distortion) 
and the complete CDWBA framework are compared with the predictions adopting the EMA. 

The gross features which were pointed out in the previous sections again show up. For 
the two proton kinetic energies, electron distortion shifts the reduced cross section towards 
higher missing momenta. However this shift is less pronounced in the CDWBA calculation 
than in the EMA approach. In Table |V| we list the missing momenta corresponding with 
the first maxima in the reduced cross section for knockout from the 2pl/2 orbit for the 
different steps in the formalism. We remark a general behaviour for the two proton kinetic 
energies. Including FSI effects which is equivalent with going from a PWIA to a DWBA 
approach, a shift towards lower p m is noticed. This shift, opposite to the shift due to 
electron distortion, can be easily explained on the basis of an effective proton momentum. 
The ejected proton feels an attractive potential (real part of the optical potential) which 
causes the detected proton to have a smaller asymptotic momentum p p than the momentum 
Pp of the initially struck proton. Table [V] also shows that this shift towards lower p m is 
increasing with decreasing proton energy. The latter is easily explained as the real part of the 
optical potential induces a shift in the average measured proton momentum approximately 
given by fill 
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^( 1 + f <v> )^ 



(38) 



where < V > is the average value of the real part of the optical potential over the interaction 
region. 

The inclusion of electron distortion effects in the model shifts the reduced cross section 
towards higher p m . The shift obtained from the EMA is larger for T p = 100 MeV than for 
T p = 70 MeV since for the latter the reduced cross section at the peak position corresponds 
with a smaller momentum transfer. The complete CDWBA calculation produces more than 
just a shift towards higher p m . The focusing of the electron beam in the vicinity of the 
target nucleus strongly modifies the maxima and minima of the reduced cross section with 
respect to the DWBA results. Since the extracted spectroscopic factors are sensitive to the 
behaviour of the reduced cross section at the peaks, an accurate prediction of this focusing 
effect is extremely important for an accurate deduction of these quantities. 

In Fig. |6] the CDWBA results are confronted with the data. The different curves are 
multiplied with the spectroscopic factors that are determined from the T p = 70 MeV data 
(Table |TV| ). Firstly, it is clear that for T p = 70 MeV the calculated cross sections are in very 
good agreement with the data. On the other hand, for the T p = 100 MeV data the shape of 
the measured reduced cross sections is not well reproduced by the CDWBA calculations and, 
as such, the extracted spectroscopic factors can not be considered as reliable. This conclusion 
agrees with the findings of den Herder in ref. [ 20fl . In ref. |20|| it was shown that a slight 



reduction of the depth of the central imaginary part of the optical potential resulted in a 
much better agreement with the data for T p = 100 MeV. With this modified optical potential 
an equally good fit of the elastic proton scattering data was obtained. This indicates that 
low energy (p,p') reactions are rather insensitive to the depth of the imaginary part in 
the nuclear interior. The (e,e'p) results, however, are sensitive to this part of the optical 
potential. The second maximum in the 2p reduced cross section reflects the behaviour of the 
2p single-particle wave function in the nuclear interior. Since the overlap is taken with the 
continuum wave functions, the second maximum is sensitive to the shape of the continuum 
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wave function in the nuclear interior. Given the uncertainties in the optical potential, this 
behaviour is not very accurately determined. As the 1/5/2 single-particle wave function 
is more surface peaked, this also explains why the 1/5/2 reduced cross section is not that 
sensitive to the depth of the imaginary part of the optical potential. Clearly, the sensitivity 
of the reduced cross section to the parameterization of the optical potential is a general 
weakness of CDWBA models but does not affect the general conclusions with respect to the 
role of electron distortion on the reduced cross section. 

The spectroscopic factors extracted from our 90 Zr(e,e'p) calculation are systematically 
larger than the values obtained by den Herder [20| . This deviation can be partly attributed 
to the fact that in ref. |2(J a different optical potential is considered. Furthermore, the 
analysis performed by den Herder accounts for electron distortion effects in an approximate 
way, thus overestimating the focusing effect of the electron beam. 

The electro-induced one-proton knockout reaction from 208 Pb is the ultimate testing 
case to study electron distortion effects. The 82 protons in 208 Pb generate a strong Coulomb 
potential felt by the initial and final electron. We have calculated the 208 Pb(e,e'p) reduced 
cross sections for proton emission from the 3sl/2, 2d3/2, 2d5/2, lg7/2 and l/ill/2 shells. 

As the effect of Coulomb distortions increases with proton number Z, we consider 208 Pb 
the ideal target nucleus to illustrate the numerical accuracy of our technique. The conver- 
gence rate for the electron partial waves is illustrated in Fig. |7|. Convergence is reached for 
I = 50 and the code is verified to produce gradually converging results, which is not evident, 
considering the large number of partial waves that has to be considered. The convergence 
tests were performed with electron wave functions of the spherical Bessel type. Accordingly, 
when convergence is reached the resulting cross section should coincide with the one ob- 
tained in a DWBA approach, provided that similar model assumptions with respect to the 
bound state wave functions and the FSI are adopted. In the insert of Fig. [7| it is verified 
that our CDWBA code bears this thorough test. 

Fig. R] shows a comparison of the CDWBA reduced cross sections with the DWBA and 



EMA results. All curves are multiplied with the spectroscopic factors as listed in Table [Vl 
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The spectroscopic factors are derived from a best fit of the CDWBA curve to the data. The 
data are well reproduced in the CDWBA, especially for the positive p m side. We remark that 
electron distortion effects considerably improve on the agreement with the data. Besides a 
shift towards higher p m , the minima and the maxima of the DWBA reduced cross section 
are strongly modified when including electron distortion effects in the model. 

In Table |VI1| some characteristics of the first and second maximum of the 3sl/2 reduced 



cross section are listed. The numbers clearly demonstrate that the final-state interaction 
of the ejected nucleon with the residual nucleus causes a small shift of the reduced cross 
section towards higher p m . Electron distortion effects also show up in a shift towards higher 
p m . The EMA again overestimates this feature compared with the complete distorted wave 
calculation. Moreover, in the CDWBA model, the shift in p m related to electron distortion 
is more pronounced for the first than for the second peak. This can be easily understood by 
considering that the two peaks in the 3sl/2 reduced cross section correspond with different 
values for the momentum transfer q. The focusing of the electron beam onto the target 
nucleus is reflected in an enhancement of the cross section at the peak positions with respect 
to the EMA cross section. The relativistic calculation by Udias et al. |T2| predicts a relative 
enhancement which is somewhat larger than our estimate. Nevertheless, it is clear that 
the two complete distorted wave calculations do not reproduce the strong focusing effect of 
electron distortion as observed with the DWEEPY code |p8| , ^9| . The CDWBA model of the 
Pavia group predicts an enhancement for the first peak in the 3sl/2 reduced cross section of 
about 20% due to the focusing of the electron beam onto the nucleus. This model accounts 
for electron distortion effects up to second order in the high-energy approximation. It has 
to be stressed that apart from the treatment of electron distortion effects, the Pavia and our 
model are very similar. Accordingly, the procedure of treating electron distortion effects is 
the only plausible explanation for the considerably different spectroscopic factors extracted 
with the two models. 

Summarizing, the effect of Coulomb distortion in parallel kinematics can be understood 
in terms of a shift of the reduced cross section towards higher p m and a small enhancement 
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of the reduced cross section due to the focusing of the electron beam onto the nucleus. It is 
also demonstrated that electron distortion effects become more important for heavier nuclei 
and need to be treated in a complete distorted wave calculation in order to extract reliable 
spectroscopic factors. 

B. Constant q — u kinematics 

We find that for constant q — u kinematics electron distortion effects act in a different 
way than for parallel kinematics. This is illustrated in Fig. || where we investigate electro- 
induced proton knockout from the 3sl/2 shell in 208 Pb for quasielastic kinematics (p p ~ q). 
In this calculation, the electron energy is the same as the one for the parallel kinematics case 
considered in the previous section. The EMA no longer causes a shift of the reduced cross 
section towards higher p m but now results in a strong quenching of the reduced cross section 
around p m = 0. The peak of the reduced cross section at p m = in Fig. |9] coincides with 
the reduced cross section for parallel kinematics at p m = displayed in Fig. |[ Moreover, 
the reduced cross section in parallel kinematics around p m = shows a strong dependence 
on q, i.e., the slightest modification of q considerably affects the value for the reduced cross 
section. Consequently, a different prescription of q than the one extracted from eq. fl2"5[) 
can cause a considerable modification of the EMA reduced cross section at p m = for 
constant q — u kinematics. 

In going from the EMA to the CDWBA, a strong enhancement of the 3sl/2 reduced 
cross section around p m = MeV/c is observed. From the previous considerations, this 
can be ascribed to a smaller g eff -value than the one adopted in the EMA approach and the 
focusing effect of the electron beam onto the nucleus. For constant q — uo kinematics it is 
therefore more difficult to disentangle the different contributions from electron distortion as 
the momentum transfer shift and the focusing effect both might cause either an enhancement 
or quenching of the reduced cross section. 

In order to give a complete picture of the role of Coulomb distortion on the exclusive 
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(e,e'p) reaction, the 208 Pb reduced cross section for constant q — u kinematics is also studied 
for non-quasielastic kinematics. Two different kinematics are considered, one in the low- 
(g > p p ) and one in the high-energy (q < p p ) side of the quasielastic peak. The results are 



plotted in Fig. 10. 



The reduced cross sections show a specific pattern regarding electron distortion effects. 
The DWBA results are shifted towards higher (lower) missing momenta for q < p p (q > p p ). 
This feature can be easily explained within the EMA. We stress that this EMA shift is 
only a first estimate of the role of electron distortion on the exclusive (e,e'p) reduced cross 



section for non-quasielastic kinematics. It is clear from Fig. 10 that a complete distorted 



wave calculation is required to completely account for electron distortion effects. 

V. CONCLUSIONS 

We have analyzed the quasielastic (e,e'p) reaction from various target nuclei. A technique 
is presented to deal with Coulomb distortion effects in an exact manner keeping the compu- 
tational time within reasonable limits. We presented results for reduced cross sections as a 
function of missing momentum, corresponding to proton knockout from the outermost shells 
in 16 0, 40 Ca, 90 Zr and 208 Pb considering two different types of kinematical arrangements. 

For parallel kinematics, Coulomb electron distortion causes a shift of the cross section 
towards higher missing momenta. This shift can be partially reproduced by the use of an 
effective momentum transfer. Furthermore,the focusing effect of the electron beam onto the 
nucleus mainly affects the maxima and minima of the reduced cross section. It has been 
shown that in order to extract realistic spectroscopic factors an accurate determination of 
this effect is highly needed. 

The role of electron distortion in the (e,e'p) reaction for constant q — u kinematics was in- 
vestigated on the basis of the one-proton knockout reaction from the 3sl/2 shell in 208 Pb. For 
quasielastic kinematics the focusing effect and the effective momentum transfer effect caused 
by Coulomb electron distortion are reflected in an enhancement, respectively, quenching of 
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the DWBA reduced cross section around p m = 0. For non-quasielastic kinematics, electron 
distortion effects cause a shift of the DWBA reduced cross section towards higher or lower 
missing momentum depending if we probe the high- or low-energy side of the quasielastic 
peak. 

The spectroscopic factors extracted in our model agree within 20% with the correspond- 
ing spectroscopic factors derived from an analysis applying the CDWBA code of the Pavia 



group. As can be seen from Fig. [11], we mostly obtain larger values. This can be partially 
attributed to the way of treating electron distortion effects. The approximate treatment of 
electron distortion by the Pavia group JTIJ tends to overestimate the focusing effect of the 
electron beam onto the nucleus and produces in this way smaller spectroscopic factors than 
our values. 
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APPENDIX A: FEYNMAN AMPLITUDE IN THE CDWBA 

In this appendix we work out the CDWBA feynman amplitude (0) for the electro- 
induced one-nucleon knockout process. Combining eqs. fl32|) and (|35|) with the expressions 
for the electron charge-current four- vector (f|), the distorted electron wave functions (|2l|) and 
the distorted nuclear wave function ( pl|) and applying some basic properties of the Wigner 
T> matrices, the summation (RBI) can be rewritten as 
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For the derivation of this expression we have chosen the reference frame (x, y, z) and the 
electron and nucleon scattering angles according to the definitions fixed in Fig. [l](b). 

The electron part £ reads as 
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The radial integrals T^Lhn lY1 this expression are evaluated in the following way 
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where n stands for j i+^+ L i n the Coulomb and electric radial integrals and for 'i+^+ L ~ 1 j n 
the magnetic radial integrals. It can easily be verified that in case l±, l 2 satisfy l± + l 2 + L even 
(odd) than 7 1; l 2 satisfy li+l 2 +L even (odd). Moreover, RL,j!,j 2 ( e J e '? <?) can be evaluated with 
either of the two choices (l\, l 2 ) or (Ji,J 2 ) since for both sets the same result is obtained. The 
initial and final radial electron partial waves occur in the expressions for the radial integrals 
Rinhhh given by 
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The following reduced matrixelements are useful in the evaluation of these integrals 
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The nuclear part M in eq. ( |A^ ) is written in terms of the reduced matrixelements £ cou b 
C e \ and £ mas defined as follows: 
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We get 
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At this stage we described the (e,e'N) process in its most general form. All approxima- 
tions with respect to the photoabsorption mechanism and the final-state interaction (FSI) 
of the ejected nucleon with the nucleus are contained in the matrixelements C cou i(C; qooL), 
C c \(C;qujL) and £ mag (C; quiL). Moreover, besides the fact that we consider the ultrarela- 
tivistic limit, electron distortion effects are accounted for exactly. For the one-body nuclear 
current operator of the impulse approximation (|26|), these reduced matrixelements are eval- 



uated in ref. [30 . 
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TABLES 





e (MeV) 


uj (MeV) 


p m (MeV/c) 


16 0(e,e'p) M 


455.8 


115 


-177 ... 265 


40 Ca(e,e'p) W 


460 


114 


-225 ... 285 


90 Zr(e,e'p) W 


346.5 


81 


27 . . . 168 


90 Zr(e,e'p) ( c ) 


350.7 


114 


62 . . . 298 


208 Pb(e,e'p) W 


412.3 


113 


-50 . . . 300 



to ref. |3; W refs. ||,|T| M refs. J20],|2T|; ^ ref. 

TABLE I. Kinematical conditions for the considered reactions. 





£* (MeV) 


CDWBA (this work) 


CDWBA^ 


lPl/2 
1P3/2 



6.3 


0.66 
0.54 


0.64 

0.51 



( a ) ref. fli~7 



TABLE II. Spectroscopic factors for the 0(e,e'p) N reaction 





^ (MeV) 


CDWBA (this work) 


CDWBA^ 


rel. CDWBA* 6 ) 


rel. CDWBA^ 


1^3/2 

2si/ 2 



2.522 


0.60 
0.48 


0.65 
0.51 


0.80 

0.75 


0.76 (0.60) 
0.51(0.44) 



to ref. 0; W ref. [0; W ref. ^2\ 

TABLE III. Spectroscopic factors for the Ca(e,e'p) 39 K reaction. The spectroscopic factors 

between brackets are obtained with the ccl nuclear current operator instead of the cc2 operator 

usually adopted in the relativistic calculations. 
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£* (MeV) 


CDWBA (this work) 
(T p = 70 MeV) 


CDWBA (this work) 
(T p = 100 MeV) 


CDWBA^ 


2pi/2 





0.42 


0.31 


0.34 


2p 3 /2 


1.507 


0.51 


0.36 


0.44 


1/5/2 


1.745 


0.52 


0.44 


0.33 



(a) ref _ 



TABLE IV. Spectroscopic factors for the Zr(e,e'p) Y reaction 



„'r^8EH 





p m (MeV/c) 


p m (MeV/c) 




T p = 70 MeV 


T p = 100 MeV 


PWIA 


61 


61 


DWBA 


50 


56 


EMA 


66 


77 


CDWBA 


61 


66 



TABLE V. The missing momentum corresponding with the first peak in the 2pl/2 reduced 
cross section for the different approaches. 
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E x (MeV) 


CDWBA (this work) 


CDWBA^ 


rel. CDWBA W 


rel. CDWBAW 


3si/ 2 





0.51 


0.40 


0.71 


0.70 (0.65) 


2^3/2 


0.35 


0.54 


0.46 




0.73 (0.66) 


2^5/2 


1.67 


0.41 


0.39 




0.60 


1^11/2 


1.35 


0.43 


0.42 




0.64 


lflY/2 


3.47 


0.21 


0.19 




0.30 



( a ) ref. ||; W ref. [0; M ref. @^ 

TABLE VI. Spectroscopic factors for the 208 Pb(e,e'p) 207 Tl reaction. 





Pm (MeV/c) 




focusing effect 






first peak 


second peak 


first peak 


second peak 


PWIA 







2.04 




DWBA 


3 


195 


1.00 (1.00 ; 1.00) 


1.00 (1.00 ; 1.00) 


EMA 


32 


209 


0.96 


1.04 


CDWBA 


24 


200 


0.99 (1.08 ; 1.21) 


1.08 (1.14 ; -) 



TABLE VII. The missing momentum and the value of the reduced cross section relative to the 
DWBA result corresponding with the first and second peak of the 3sl/2 reduced cross section. 
The corresponding values obtained by Udias et al. [12] and Giusti et al. [28,^9| are listed between 
brackets. 
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FIGURES 




(a) 




(b) 



FIG. 1. Kinematics for the (e,e'N) reaction in the CDWBA. 
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Electron part 



Nuclear part 









n%^ 2 (e,e';q) 1/2 < j, < j liTaax 
Tli jlj2 (e,e';q) < L < L max 

nl Ljij 2 ( e > e '; q) ?min < g < g ma x 



^coui(C;g^L) 
£ e i(C; gcuL) 
£ma,g{C:, quiL) 


S -^ S ^max 
^min _; q _: Qmax 



integral over q 


sum over L and j\ 


U 


™CDWBA 



FIG. 2. Schematic representation of the CDWBA approach for the exclusive (e,e'p) cross section. 
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100 200 

P m (MeV/c) 

FIG. 3. Comparison of the DWBA (dashed line) and CDWBA (solid line) results for pro- 
ton knockout from 16 for parallel kinematics. The curves are multiplied with the appropriate 
spectroscopic factors from Table [D]. The data are from ref. [17]. 
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FIG. 4. Comparison of the DWBA (dashed line) and CDWBA (solid line) results for proton 
knockout from 40 Ca under parallel kinematics. The curves are multiplied with the appropriate 



spectroscopic factors (see Table HI). The data are from ref. M 
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T p = 100 MeV 



2p 1/2 (xlO z ) \ 




10 
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100 200 
P m (MeV/c) 



FIG. 5. Effect of Coulomb distortion on the reduced cross sections for proton knockout from 
the three valence shells in 90 Zr at T p = 70 MeV and T v = 100 MeV. The dashed line stands for the 
DWBA result, the dotted line the EMA result and the solid line the complete CDWBA calculation. 
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FIG. 6. Comparison of the CDWBA calculation with the NIKHEF data for proton knockout 
from 90 Zr ( f20|) (solid line: T p = 70 MeV; dotted-dashed line: T p = 100 MeV). The curves are 
multiplied with the spectroscopic factors derived for T p = 70 MeV (Table [TV|). 
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FIG. 7. Convergence check of the CDWBA 208 Pb(e,e'p) calculation in parallel kinematics. The 
electron wave functions are described by spherical Bessel functions. For the dotted, dot-dashed 
and dashed line electron partial waves up to I = 30,40 and 50 are considered. In the insert the 
DWBA calculation (solid line) is compared with the CDWBA calculation when convergence is 
reached (dashed line). 
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FIG. 8. The reduced cross sections for electro-induced one-proton knockout from the valence 
shells in 208 Pb for parallel kinematics. The dashed, dotted and solid curve give the DWBA, EMA 
and CDWBA results. The calculations are compared with the data from ref. |J and are multiplied 



with the spectroscopic factors from Table. VI. 
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FIG. 9. The reduced cross section for proton knockout from the 3sl/2 shell in 208 Pb for constant 
q - to kinematics (e = 412.3 MeV, q = AAA MeV/c, to = 113 MeV). The dotted-dashed, dashed, 
dotted and solid line represent the PWIA, DWBA, EMA and CDWBA results. The curves are not 
multiplied with a spectroscopic factor. 
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-300 -200 -100 100 200 300 

p m (MeV/c) 

FIG. 10. The reduced cross section for proton knockout from the 3sl/2 shell in 208 Pb for 
constant q — to kinematics. Upper figure: e = 412.3 MeV, q = 350 MeV/c, u> = 113 MeV; bottom 
figure: e = 412.3 MeV, q = 600 MeV/c, u = 113 MeV. The dashed and solid line represent the 
DWBA and CDWBA results. The curves are not multiplied with a spectroscopic factor. 
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10 

mass number A 



FIG. 11. Spectroscopic factors derived from the A(e,e'p) reaction to the groundstate of the 
residual nucleus. The black squares give the results within the presented model, whereas the open 



squares are the values obtained with the DWEEPY code [11| which incorporates electron distortion 
effects in an approximate manner. 
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